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We derive an approximate local density functional for the exchange-correlation energy to be used 
in density-functional calculations of two-dimensional systems. In the derivation we employ the Colle- 
Salvetti wave function within the scheme of Salvetti and Montagnani [Phys. Rev. A 63, 052109 
(2001)] to satisfy the sum rule for the exchange-correlation hole. We apply the functional for the 
two-dimensional homogeneous electron gas as well as to a set of quantum dots and find a very good 
agreement with exact reference data. 
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Development in modern technology has enabled the 
fabrication of nanoscale electronic devices with a large va- 
riety of lovir-dimensional systems. Two-dimensional (2D) 
quantum dots (QDs) are particularly interesting exam- 
ples due to the tunability in their size and shape, and 
number of confined electrons P, . From the theoretical 
point of view, QDs constitute an ideal platform to study 
the many-particle problem, electronic correlations, and 
the role of the dimensionality. 

In density-functional theory Q (DFT) particle-particle 
interactions beyond the classical (Hartree) term are cap- 
tured through the exchange-correlation (xc) functional, 
which is approximated in practice. The development of 
xc functionals of varying portions of simplicity and ac- 
curacy have a long and successful history [J] . The Colle- 
Salvetti (CS) scheme H Q and its variants Vn have had 
an important role in the development, especially in terms 
of the electronic correlation. However, these efforts have 
focused almost solely on three dimensions (3D), apart 
from orbital functionals where the aspect of dimension- 
ality is inbuilt through the Kohn-Sham orbitals. Only 
very recently, several local @, Q and semi-local [Tol - [l5| 
functionals have been developed in 2D, and in many test 
cases involving, e.g., different QDs, they have outper- 
formed the commonly used 2D local-density approxima- 
tion based on the exact exchange and correlation of the 
homogeneous 2D electron gas [3; [13 (2DEG). 

In Ref. Q a 2D local density functional for the corre- 
lation energy was derived using the CS framework with 
a Gaussian summation for the pair density [isj . Despite 
the good performance of this functional for the correla- 
tion, a compatible approximation for the exchange energy 
is needed in view of, e.g., total-energy calculations. In 
fact, a combination with the 20 generalized-gradient ap- 
proximation for the exchange 113, leads to a reasonable 
accuracy in the total energy |19| . However, this com- 
bined functional is still semi- local, i.e., it depends on the 



density gradients, which reduces the numerical efficiency. 

In this work we employ the CS framework to derive 
a 2D local functional for the xc energy, so that both the 
exchange and correlation are treated in the same footing. 
In the derivation we follow the 3D scheme of Salvetti and 
Montagnani [23| for the second-order density matrix to 
satisfy the sum rule of the xc hole, which is used to obtain 
a local density functional of a simple polynomial form. 
We optimize two remaining parameters of the functional 
by fitting against exact results for six-electron QDs. The 
obtained parameters show universality in the sense that 
a good accuracy and consistency is found when testing 
the functional for QDs with a varying electron number 
as weU as for the 2DEG. 

The electron-electron interaction energy can be for- 
mally expressed (in Hartree atomic units) as 
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is the diagonal element of the spinless second-order den- 
sity matrix describing the distribution density of elec- 
tron pairs. Here / dN denotes the spatial integration 
and spin summation over the A^th spatial and spin co- 
ordinates rjvCTAr, and ^'(1, 2, . . . , A^) is the ground-state 
many-body wave function. The clement p2(vi,T2) satis- 
fies the normalization 
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= / dri / dr2/92(ri,r2), 
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and it is related with the electron density, i.e., the diag- 
onal term of the first-order density matrix through 
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P(ri) = j;^^ J dr2 P2(ri,r2). 
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Next, introducing a symmetric function accounting for 
all nonclassical effects called the pair correlation function 
/i(ri,r2) suggests us to write [22] 

P2(ri,r2) = ip(ri)p(r2)[l + /i(ri,r2)]. (5) 

The important sum rule for the xc hole can be expressed 
in terms of the pair correlation function as 



j dv2p[r2)h{ri,Y2) = -1. 



(6) 



The total interaction energy in Eq. ([T]) can be split in the 
classical Coulomb (Hartree) term associated with p(r) 
and the nonclassical (indirect) part associated with the 
xc energy, 
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It should be noted that, compared with the conven- 
tional DFT formalism, E^c as defined here neglects the 
kinetic-energy contribution to the correlation energy. 
The Hartree term can be computed in a straightforward 
fashion, but the integration in the xc part is nontrivial 
due to the pair correlation function ft,(ri,r2). The key 
point in the present work is to obtain an approximation 
for /i(ri,r2) satisfying the sum rule in Eq. Before 
proceeding with that, we will briefly introduce the CS 
approach which is relevant for the derivation. 

The CS scheme starts with the following ansatz for the 
many-body wave function Q 

*(ricri, riyaN) = *HF(ricri, tnctn) J|[l-(/3(rj, r^)], 

i>j 

(8) 

where HF refers to the single-determinant Hartrec-Fock 
wave function, and 
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with r = |ri — r2|. Above, 7 and a are optimizable pa- 
rameters (7 with dimension of p~^), and $ is a monotonic 
function varying between zero and one. The differences 
from the original CS scheme are obvious; most impor- 
tantly, /3 is now a non-local functional of the density. 
We may now search for the pair correlation function 
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where / is assumed to be a simple polynomial of the form 

/($) = ao*" + ai$""' + ... (14) 

The nominator in the expression for /i(ri,r2) is similar 
to the CS functional 0, Q , whereas the denominator is 
chosen such that to the sum rule in Eq. (jG]) is satisfied. 
Substituting Eq. dTS]) into Eq. ^ yields 

j dv2 p{Y2)h{ri,r2) = J drp{ri + r)h{ri,ri + r) 
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-2$3e'2^*r^(l + r) + $2e-/3'-^ 
X fe-'^'^' -H2(l + r)l --2$e-''''' I = -1 (15) 



This expression involves integrals which can be computed 
by using mean value theorem and the regularit y o f the 
functions. By following the procedure of Ref. [201, 
obtain 



^(ri,r2) - [1-$(R)(1 + Cr)] cxp[-/32(R)r2] (9) j dr ^(ri, n + r) e-^^'-^^'-^+'-^'-'r" 



describes the correlated part of the wave function writ- 
ten in center-of-mass, R = (ri + r2)/2, and relative, 
r = ri — r2 coordinates. The parameter C, comes from 
the cusp conditions, and the quantities $ and /3 act as 
correlation factors. In Refs. [31 and @ dealing with 
3D and 2D systems, respectively, /3 was introduced as a 
local factor for the correlation length, /3(R) = gp^/^(R), 
where D is the dimension, q is a fitting parameter, and 
p(R) is the electron density. The CS approach assumes 
that the first and second-order density matrices can be 
written as pi(ri,r2) = pf^(ri,r2) and p^^(ri,r2) = 
p^^(ri,r2) [1 - ^{Yi,V2)f, respectively ^\. 

To approximate /i(ri , r?), we extend the strategy by 
Salvetti and Montagnani [201 2D by introducing the 
correlation factors in the following way: 



27rg(ri) j dre"'"'" r"+i. (16) 



Utilizing this approximate integration, which becomes 
more accurate as 6 becomes large, leads to 



c?r2p(r2)/i(ri,r2) 
27rp(ri) 



/3(ri,r2)=7p'/'(ri)p^/'(r2), 



(10) 



$4(io + 2ii -f' i2) - $^(2io + 2n) 
+ $2(^0 + 2jo + 2ji) - 2$jo] = -1, (17) 
where we define in and j„ as 

i„ = ^2-"-2/3-" / dx e-^'a;"+i (18) 
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and 



Jn = VF" 



(19) 



Using the definition of /3 and calculating the integrals 
leads to the final result for the polynomial function, 
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where the coefficients are given by 
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Computation of the integral in Eq. ([7]) is performed by 
a similar procedure, 
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Using Eq. ([20)) and definition of /3 we obtain the final 
result for the xc energy in 2D, 



with 



_ Frr^ f bo^^ + 61$^ + 62<I> - 1 
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(24) 
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The remaining task is to find a reasonable pair of val- 
ues for 7 and a which determine q{p) through Eqs. pop . 
(|TT]) . ((2T|) . and ((23)) . Here we choose to fit these pa- 
rameters to reproduce the xc energies of parabolic QDs 
with N = 6, which is the smallest closed-shell sys- 
tem beyond the simplest N = 2 case. In the exter- 
nal potential v^xti^) = w^r^/2 we use the confinement 
strengths = 1/4 and 1/16 for which numerically ex- 
act configuration-interaction (CI) data is available [23j . 
These confinements have a rather wide range with respect 



TABLE I: Exchange-correlation energies for parabolic quan- 
tum dots. The optimal parameters 7 = 1.12 and a = 0.45 
have been used for the calculation of E^"'^"^. The last row 
contains the mean percentage error. 
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0.59871" 
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1/6 
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1/16 
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0.2789 


0.2820 
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1/1.89^ 


2.156+ 


2.138 


2.137 
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1/4 


2.014+ 


2.008 


2.011 
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1/16 


0.9265+ 


0.9309 


0.9429 


12 


1/1.89^ 


4.708* 


4.716 


4.701 


A 






1.86% 


2.19% 



* Total energy from the analytic solution in Ref. 
+ Total energy from the CI data in Ref. [231. 
+ Total energy from the QMC data in Ref. 



to the relative weight of the xc effects, and, moreover, the 
chosen values are realistic regarding the modeling of real 
QD devices [H, 0. The reference xc energy is obtained 
from 
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where fi'^of'^* is the reference total energy, Ef^^ is the 
total energy from the exact-exchange (EXX) calculation 
performed here within the Krieger-Li-Iafrate approxima- 
tion [iJl and using the octopus code [l^l, and E^^^ is 
the exchange energy. The best fit with i?^°+ is obtained 
with parameter values 7 = 1.12 and a ~ 0.45. 

Next we test if the chosen parameter values yield 
reasonable and consistent results for different 2D sys- 
tems. This is naturally a desired property in any den- 
sity functional in order to be a predictive approxima- 
tion. First we consider parabolic QDs with = 2 ... 12 
and bj = 1/16 ... 1. The results are summarized in Ta- 
blelH The reference xc energies have been calculated from 
Eg. ([26l) using the total-energy data from analytic [26l |. 
CI [23|, and quantum Monte Carlo [l^] (QMC) calcula- 
tions as indicated in the table. Overall, we find a very 
good performance of our functional, the mean error be- 
ing 1.86%, which is smaller than that of the 2D-LDA 
(2.19%). Although also the LDA is this accurate for the 
xc energy, it should be noted that both the exchange 
and correlation parts, respectively, have significant er- 
rors (sec, e.g., Refs. [l^l and [111), and the good overall 
performance follows from the well-known error cancella- 
tion. It is also noteworthy that the 12-electron case is 
very accurate, and it can be expected that the accuracy 
is preserved for larger when TV is increased further. 

Table [I] raises a natural question whether the good per- 
formance simply follows from the fact that 7 and a were 
fitted to a similar system with N = 6. Therefore, in 
Figs. [Ua) and (b) we examine the "extreme" cases of 
TableUwith TV 2 (w = 1) and = 12 (w = 1.89^), re- 
spectively. The figures show the absolute relative errors 
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FIG. 1: (color online). Absolute relative error for parabolic 
quantum dots, (a) TV = 2, a; = 1 and (b) N = 12, u = 1/1.89^ 
with respect to parameters 7 and a. The crosses mark the 
values chosen from the fit to the N = 6 case. 



as functions of both 7 and a, so that here the parameter 
values have been left undetermined for both two cases. 
The white crosses show the chosen values 7 = 1.12 and 
a = 0.45 based on the = 6 data. In both cases, the 
crosses match very well with the optimal regime where 
the relative error is smallest for iV = 2 and N ~ 12. 
Hence, Fig. [T] confirms that, at least for parabolic QDs, 
the functional is consistent. The figure also demonstrates 
the strong correlation between the two parameters as well 
as the uniqueness between them - for each 7 (a) there is 
only one compatible a (7). 

Finally we consider the 2DEG corresponding to the im- 
portant limit of an infinite electron number. Figure [H^a) 
shows the comparison of the xc energy per particle with 
respect to the exact 2DEG result. Here we have used the 
same parameter values 7 — 1.12 and a = 0.45 as before. 
We find an excellent agreement through a wide range of 
the density parameter rg = (7rp)~^/^. In the realistic 
density range the relative error is within a few percent 
(see the inset). In Fig. ^h) we show the function q{p) 
of our functional (dashed line) in comparison with the 
optimal values to reproduce the exact xc energy of the 
2DEG. Overall, we find good consistency in the results 
at varying . More importantly, regarding the values for 
7 and a the present functional is also consistent in the 
comparison between 2DEG and the QDs above. 

To summarize, we have used the Colic- Sal vetti scheme. 
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FIG. 2: (color online), (a) Exchange-correlation energy per 
particle for the two-dimensional homogeneous electron gas 
obtained using our functional with the chosen parameters 
7 = 1.12 and a = 0.45 (dashed line) in comparison with 
the exact result (solid line). The inset shows the relative er- 
ror, (b) Result for the function q{p) obtained using our func- 
tional (dashed line) in comparison with the optimal values 
required to reproduce the exact exchange-correlation energy 
of the two-dimensional homogeneous electron gas. 



and in particular its recent generalizations to derive an 
approximate local density functional for the exchange- 
correlation energy of electrons in two dimensions. The 
functional has a simple polynomial form and it fulfills 
the sum-rule constraint of the exchange-correlation hole. 
We have fitted the remaining free parameters against ex- 
act results for six-electron quantum dots and found an 
excellent consistency in the results for a set of quan- 
tum dots with varying electron number and varying rel- 
ative proportion of the exchange-correlation energy. The 
functional is precise also for the two-dimensional homo- 
geneous electron gas with the same fixed parameters. 
Therefore, we may expect the functional to have predic- 
tive power in density-functional calculations for various 
two-dimensional electron systems. In this respect, gen- 
eralization to spin-polarized systems would be the most 
important future extension of the method. 
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